Spectral function of transverse spin fluctuations in an antiferromagnet 
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The spectral function of transverse spin fluctuations, including the contributions from both the 
single-particle and the collective (magnon) excitations in an antiferromagnet, is evaluated for the 
Hubbard model with NN and NNN hoppings in the full t/-range from weak coupling to strong 
coupling. For the NN hopping model, the magnon excitations are dominant for U > 2.5 (d — 2), so 
that an effective spin description of the AF state holds down to a surprisingly low U value. In the 
weak coupling limit the spectral function is suppressed at low energies, as if due to an effective gap. 
NNN hopping t' leads to magnon softening and also a significant increase in the low-energy spectral 
function due to single-particle excitations. Evolution of the magnon spectrum with t' is studied in 
the strong coupling limit, and the quantum spin-fluctuation correction to sublattice magnetization 
in d — 2 and the Neel temperature in d = 3 are also evaluated. 
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I. INTRODUCTION 

The antiferromagnetic state of the half-filled Hubbard 
model with nearest-neighbour (NN) hopping is charac- 
terized by an energy gap, and the spectrum of trans- 
verse spin fluctuations consists of the low-lying, collec- 
tive (magnon) excitations, as well as the single-particle 
excitations across the gap. Generally, there is a clear dis- 
tinction between these two excitations which are well sep- 
arated in energy. However, recent extensions, e.g. with 
disorder, impurities, Q and next-nearest-neighbour 
(NNN) hopping, || clearly show the existence of essen- 
tially gapless antiferromagnetism, arising from a variety 
of mechanisms. Thus, with increasing disorder (on-site 
potential disorder) the two Hubbard bands progressively 
broaden, until the band gap vanishes when the disorder 
strength W » U . With \aw-U impurities on the other 
hand, the effective charge gap becomes negligible due to 
nearly localized states on the impurities. [|2| And lastly, 
when NNN hopping is included, a band asymmetry is in- 
troduced which reduces the band gap in the AF state, 
and in weak coupling there exists a region of gapless AF 
phase in the magnetic phase diagram. |J 

This gaplessness implies that the collective magnon ex- 
citations are no longer well defined, and actually merge 
with the single-particle excitations, thereby necessitat- 
ing a unified scheme for the evaluation of transverse spin 
fluctuation. While the magnon contribution to the trans- 
verse spin fluctuations was evaluated recently, and the 
sublattice magnetization and the Neel temperature were 
obtained within a renormalized spin fluctuation theory 
in the whole U/t range, Q|| in this paper we describe an 
alternate scheme for evaluating the transverse spin fluc- 
tuations which allows both the collective excitations and 
the single-particle excitations to be studied on the same 
footing. 

Single-particle excitations are especially significant be- 
cause it is precisely this part of the transverse spin spec- 
tral function which allows for a quantitative distinction 
of the antiferromagnetic state within the Hubbard model 



from that of an equivalent Heisenberg spin model with 
[/-dependent, extended-range spin couplings Jij(U), but 
possessing only magnon excitations. A study of the rel- 
ative strengths of the magnon and single-particle exci- 
tations, in terms of the integrated spectral weights in 
the whole U/t range, will therefore allow for a quantita- 
tive demarkation along the U/t axis below which single- 
particle excitations are the dominant contribution. Thus, 
while in this low-U regime the AF state is not well de- 
scribed in terms of an effective Heisenberg spin model 
possessing only magnon excitations, in the intermediate 
and strong coupling limits use of an effective Heisenberg 
model with [/-dependent spin couplings Jij(U), || and 
generally use of a spin picture ^,|) is appropriate. 

Yet another significance of the single-particle excita- 
tions is interestingly connected with the spin commuta- 
tion relation = 2S Z . The RPA-level ground- 
state expectation value ([S + , S'~])rpa, involving the dif- 
ference of transverse spin correlations evaluated in the 
AF ground state, should be identically equal to (2S z )hf> 
the local magnetization at the HF level. This is deduced 
from the fact that both RPA and HF approximations are 
of the same order (0(1)) within the systematic inverse- 
degeneracy expansion scheme Hj in powers of 1 /Af ', where 
Af is the number of orbitals per site. Therefore both RPA 
and HF become exact in the limit Af — > oo, when all cor- 
rections of order X/Af or higher vanish. The magnon 
contribution ([S + , S~]}^p^ on was indeed found to be in 
excellent agreement with the HF magnetization (2S z )-hf 
in the intermediate and strong coupling limits. J|] How- 
ever, a discrepancy was observed at small U which was 
attributed to the neglect of the single-particle excitations, 
which become relatively more important in the weak cou- 
pling limit. We will show here that indeed when the 
single-particle excitations are included this discrepancy 
is exactly removed. 

Finally, the unified scheme for incorporating both the 
single particle and collective excitations in the evaluation 
yields the complete spectrum of transverse spin excita- 
tions in the magnetic state. Calculations with realistic 
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models of magnetic solids can hence be used for com- 
parison with results of scattering studies with neutrons, 
which are direct experimental probes into the spectrum 
of magnetic excitations in solids. In this regard the ne- 
cessity of more realistic models which include NNN hop- 
ping etc. has been acknowledged recently from realistic 
band structure studies, photoemission data and neutron- 
scattering measurements of high-T c and related materi- 
als. [§-(11 1 



II. HUBBARD MODEL WITH 
NEXT-NEAREST-NEIGHBOUR HOPPING 

We consider the following Hamiltonian on a square 
lattice, with hopping terms t and t' between nearest- 
neighbour (NN) and next-nearest-neighbour (NNN) pairs 
of sites (ij) and (ik) respectively, 



XN 



NNN 



H=-t^ 4a a 3<r - l ' 4a a k<r + U^n^n^ . (1) 



{ik)c 



Extension to the simple cubic lattice is considered in the 
Appendix. In the plane- wave basis defined by ai a — 

y^"X4 e lk ' ri ai C a, the non-interacting part of the Hamil- 
tonian H = 2kcr( e k + e k) a iL a k(T, where £k and e k are 
the two free-particle energies, correponding to NN and 
NNN hopping respectively, 



ek = — 2i(cos k x + cos k y ) 
e k = — 4t' cos k x - cos k y . 



(2) 



For the NN model, the HF-level description of the 
broken-symmetry AF state, and the transverse spin fluc- 
tuations have been discussed earlier in the strong, |fj^] 
intermediate, and weak coupling |u] limits. We 

briefly discuss the extension for the NNN hopping. Since 
the NNN hopping term connects sites in the same sub- 
lattice, in the two-sublattice basis the e k term appears in 
the diagonal matrix elements of the HF Hamiltonian 



^HF(k) 



-a A -i 



s k £k 
a A + 





-o-A e k 


= 4i + 


e k a A 



(3) 



for spin a. Here 2A = mU , where m is the sublattice 
magnetization. For the NN hopping model 2A is also 
the energy gap for single-particle excitations. Since the 
e k term appears as a unit matrix, the eigenvectors of the 
HF Hamiltonian remain unchanged from the NN case, 
fl2f whereas the eigenvalues correponding to the quasi- 
particle energies are modified to, 



E, 



(±) 

k(T 



4± 



(4) 



The two signs ± refer to the two quasiparticle bands. The 
band gap is thus affected by the NNN hopping term, and 



it progressively decreases as 2A — it' in the weak coupling 
limit. 

As the eigenvectors of the HF Hamiltonian are un- 
changed, the self-consistency condition retains its form, 
and therefore the sublattice magnetization is indepen- 
dent of t 1 ', provided there is a nonzero band gap, with 
the lower band occupied and the upper band empty. We 
will restrict ourselves only to this regime where there is 
no band overlap, though the gap may vanish when the 
two bands are just touching. The fcrmionic quasiparti- 
cle amplitudes a\^ a and feker for spin a =|, J, and the two 
quasiparticle bands are given by 



Me — u kie — "kie — "kt 




(5) 



These relationships follow from the spin-sublattice and 
particle-hole symmetry in the AF state of the half-filled 
system. 

However, this situation changes when with increasing 
NNN hopping t' the two bands overlap and the band gap 
vanishes. This overlap indicates that the low-lying states 
in the upper band corresponding to double occupancy 
are lower in energy than the high-energy states in the 
lower band. This results in charge transfer and conse- 
quently partially empty and doubly occupied sites. The 
situation is therefore analogous to the addition of holes 
or electrons in the half-filled Hubbard model, with all its 
associated complications of spin bags,strings of upturned 
spins, spiral and striped phases etc. [|15[ Thus the simple 
gapless AF state of the NNN hopping model, 0] which 
is indeed a self-consistent HF solution, may actually be 
unstable due to same instabilities as the hole-doped AF. 
Numerical HF studies do show this instability, and will 
be further discussed elsewhere. |l^] 



III. THE TRANSVERSE SPIN SPECTRAL 
FUNCTION 



The spectral function for transverse spin fluctuations 
in the AF state is obtained from the imaginary part of the 
corresponding time-ordered propagator of the transverse 
spin operators S~ and Sj~ at sites i and j, x h ( ( l w ) = 

/ ^E l e M( *- t ' ) e- l< J-^- r ^(* G |r[5r(t)5+(t')]|*G)- In 
terms of the RPA expression [x h (qw)]p.pA = 
[x°(qo;)]/[l — £/x°(q w )L where x°(q w ) is the zeroth- 
order particle-hole propagator, there are two sources 
of contribution to the spectral function, A(lj) = 
J^ q Im Tr[x ^(qui)]. The contribution arising from the 
imaginary part of x°(qo;) is associated with the single- 
particle excitations across the gap (with u) > 2A for the 
NN model), whereas that due to the vanishing of the 
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denominator 1 — C/x (qw) = is associated with the col- 
lective magnon excitations (involving ui < 2 A for the NN 
model). As mentioned earlier, this latter contribution to 
the integrated spectral weight 7r _1 J dujA(ui), and there- 
fore to the transverse spin correlations (S + S~)bj>a and 
(5 _ S' + )rpa was recently studied in the context of quan- 
tum correction to sublattice magnetization and the Neel 
temperature in the whole U/t range. [Q,^) 

The evaluation of the transverse spin spectral function, 
including contributions from both the single-particle and 
collective excitations, is facilitated by expressing the 2x2 
complex matrix [x°(qw)] in terms of its eigenvalues A™(cj) 
and eigenvectors |0™(cj)), and we obtain 

AH = £ ImTr[ X -+ (q .)] = £ Iml ^J^ 



E ImTr £ 

q n=l,2 



=£*»£ 

q n=l,2 



1 



1 - UX^to) 

a» 



(6) 



If instead of taking the trace in Eq. (6), the two di- 
agonal matrix elements are considered separately, then 



we obtain the transverse spin correlations on the A(B) 
sublattice by integrating over frequency: 

(S + S-) B(A) = (S-S + ) A(B) = J ^^Im[ X - + (q^] A (B) 



y lm ( A » 



q n=l,2 



I^MIacb)- 



Here the relationship (S + S~)b(A) = S + ) A(B) be- 
tween the transverse spin correlations on opposite sub- 
lattices follows from the spin-sublattice symmetry in the 
AF state. 



IV. COMPUTATIONAL PROCEDURE AND 
RESULTS 

To begin with the zeroth order 2x2 matrix [x°(qu;)] = 
* / f^Ek[ GT ( kV )][GHk' -q,w' - w)]> is evaluated in 
the broken-symmetry AF state. In terms of the fcrmionic 
quasiparticle amplitudes and energies, [x°(qcj)] is given 

by, P| 



(7) 



[x°M] 



£ 

k 

£ 

k 



r a kTe a k- q |€ 



akTe^kTe a k-qi©&k-qie 



akTe^kTe a k- q |e 



frk-qlf 



b 2 h 2 
°kTe°k-qie 



a kT© a k-qie 
akTe^kTeak-qie^k-qie 



akTe^kT© a k-qie^k~qie 



1 



^-q ' 



b 2 



b 2 



- LO - 



IT] 



(8) 



Analytical evaluation of [x°(qw)] m the strong and inter- 
mediate coupling limits has been discussed earlier for the 
NN hopping model, [ ^2|Jl3| and is described in the Ap- 
pendix for the NNN hopping model. For arbitrary U, the 
k-sum is numerically performed to evaluate the complex 
matrix [x°(qw)], which is then diagonalized to obtain the 
two eigenvalues A™ and eigenvectors |0q). The infinites- 
imal r\ is appropriately chosen according to the fineness 
of the k and q grids. Typically, the grid sizes dk and dq 
were taken of order 0.05 and r\ ~ 0.01. Summing over 
the whole range of q values between and n in each di- 
mension then yields the spectral function A(u>) from Eq. 
(6), and the transverse spin correlations from Eq. (7). 

The integrated spectral weight is obtained by numer- 
ically integrating the spectral function A(lu) over fre- 
quency. The collective and single-particle contributions 
can be evaluated separately by integrating over the ui- 
rcgions uj < 2A and u> > 2A respectively (for the NN 
case). As the single-particle excitations have a contin- 
uum distribution in Eq. (6), the evaluation of this con- 
tribution is less computationally intensive. However, for 
collective excitations, which are a set of delta functions, 
it is necessary to have in the numerical integration pro- 
cedure fine enough u>- and g-grids, so that a sufficiently 
large number of magnon modes are picked up in the u> 



and q sums. The method described in ref. ||4|,£J is more 
efficient for evaluating the collective-excitation contribu- 
tion. 

Before discussing the results we first consider the two 
limiting cases. In the limit of vanishing interaction 
strength U — > 0, A — > 0, we have x~ + ( c L UJ ) = X°(q w ), 
and therefore the magnon contribution vanishes. Also, 
as all the quasiparticle probabilities a 2 = b 2 = 1/2 from 
Eq. (5), we have from above 



Imx°(qw) = y, 



l l 

I I 

L 4 4 



X 7T X 



W^-q " ^k + + <^k + ^k-q - ")]■ 



(9) 



In the non-interacting limit, the integrated spectral 
weight therefore yields 



Lim 

U-yQ 



du> 



]TlmTr[ X -+(qc,)] = l/2. 



(10) 



On the other hand, in the strong coupling limit U — > 
oo, the quasiparticle bands are narrowed to infinitesimal 
width, so that the imaginary part of x°(q^) strengthens 
to delta functions at the gap frequencies u> — ±U. 
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FIG. 1. The spectral function for transverse spin fluctu- FI G. 3. The near merging of the collective (w < 2A) 

ations obtained from Eq. (6), showing a distinct separation an d the single-particle (u > 2 A) contributions to the trans- 

between the low- lying, collective excitations (ui < 2 A) and verse-spin spectral function for 2 A = 1. 
the single-particle excitations (w > 2A) . 
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FIG. 2. The single-particle contribution to the trans- 
verse-spin spectral function A(uS) = ^ Im Trx~ + (fj"^), for 
gap values 2A = 1,2,3,4, showing the rapid decrease in the 
spectral weight with increasing correlation. 

Therefore, for both possible values of yp(qu), zero or in- 
finity, the single-particle contribution to the imaginary 
part of x ^(q<^) vanishes. Thus, with increasing inter- 
action strength, the spectrum of transverse spin fluctua- 
tions changes from predominantly single-particle excita- 
tions in weak coupling to predominantly magnon excita- 
tions in the strong coupling limit. 

Figures 1 through 6 depict results for the NN model. 
Fig. 1 shows a typical transverse spin spectral function 
A(w) = J2 q ^ m Trx h (q^), for a moderate correlation 
value, showing the distinct separation between the low- 
lying, collective excitations (lo < 2A) and the single- 
particle excitations (uj > 2 A). 
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FIG. 4. The integrated spectral weights for the sin- 
gle-particle (solid) and collective (dashed) excitations, show- 
ing their strong suppression respectively in the strong and 
weak coupling limits. U ~ 2.5 marks the point below which 
an effective spin description of the antiferromagnetic state 
within the Hubbard model starts to fail. 

Fig. 2 shows the single-particle contribution to the 
spectral function for moderate gap values 2A = 1, 2, 3, 4. 
The rapid decrease in the spectral weight with increasing 
correlation shows that the single-particle excitations are 
strongly suppressed in the strong correlation limit. 

With decreasing correlation, the separation between 
the low-lying collective excitations and the single-particle 
excitations progressively decreases. Fig. 3 shows the near 
merging of the collective (uj < 2 A) and the single-particle 
(uj > 2 A) excitations for 2A = 1 (U = 2.28). 

A comparison of the integrated spectral weights as a 
function of U is shown in Fig. 4, for the single-particle 
and collective excitations. While the collective excita- 
tions are seen to sharply fall-off in the weak coupling 
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limit, 
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FIG. 5. The low-energy part of the transverse-spin spec- 
tral function for very low A, showing a gap-like suppression 
of excitations at low energies. 



FIG. 7. Effects of the NNN hopping t' on the transverse 
spin spectral function, showing the shift in the magnon spec- 
trum towards lower energy with increasing t' . 
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FIG. 6. The drastic suppression in the contribution of 
the low-energy, collective excitations in going from U = 3.29 
(A = 1.0) to U= 1.19 (A = 0.1). 

the single-particle excitations are likewise suppressed in 
the strong coupling limit. For U > 2.5, the collective 
(magnon) excitations are dominant, and thus an effective 
spin description of the antiferromagnetic state is appro- 
priate down to a surprisingly low U value. 

We now examine the role of the single-particle excita- 
tions in the sum rule following from the spin commu- 
tation relation [S' + ,5' _ ] = 2S Z . As discussed in the 
Introduction, the RPA-level, AF ground-state expecta- 
tion value ([S^, S~])rpa, including both magnon and 
single-particle contributions, should be identically equal 
to (2S()hf. The single-particle contribution to this ex- 
pectation value, obtained from the transverse spin corre- 
lations (S + S~) and (S~S + ) evaluated from Eq. (7), are 
shown in Table I for several U values. Also shown are the 
magnon contributions which were obtained earlier. [QJ^] 
The total of the magnon and single-particle contributions 
is indeed in close agreement with the HF magnetization 
(25?>hf. 



The sharp fall-off of the gapless magnon contribution 
in the weak coupling limit essentially leaves only the 
single-particle contribution. As these excitations have a 
minimum-energy threshold of 2A, the spectrum of trans- 
verse spin fluctuations essentially acquires a pseudo gap 
at low energies. This is clearly seen in Fig. 5, showing a 
gap-like suppression in the density of excitations at low 
energies. 

Fig. 6 shows a comparison of the spectral functions 
for A = 0.1 and A = 1.0. The drastic suppression in the 
contribution of the low-energy, collective excitations in 
going from U = 3.29 (A = 1.0) to U = 1.19 (A = 0.1) is 
quite remarkable. 

Fig. 7 shows the effects of the NNN hopping on the 
spectral function. With increasing NNN hopping t' , the 
magnon spectrum shifts towards lower energy, which re- 
flects the magnon softening. There is also a substantial 
increase in the low-energy spectral function due to the 
reduction in the energy gap with NNN hopping. For 
U = 3.29, we have A = 1, and the energy gap 2A — At' 
just vanishes when t' = 0.5. It is therefore clear that the 
strong enhancement in the low-energy spectral function 
for t' = 0.5 is due to the gapless single-particle excita- 
tions. 

TABLE I. The RPA-level expectation values of the spin 
commutator [Sf ,ST~], evaluated in the AF state from the 
transverse spin propagator, with contributions from the col- 
lective (magnon) excitations, and the single-particle excita- 
tions. The total is compared with the HF magnetization 
(2S ! ? )hf- 



A 


U 


(K,s-]) 


mag {[S^, Sf ]) sp 




(2S?)hf 


0.6 


2.499 


.315 


.172 


.487 


.480 


1.0 


3.292 


.553 


.064 


.617 


.607 


2.0 


5.127 


.776 


.009 


.785 


.781 


3.0 


6.948 


.865 






.863 
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This magnon softening has been analytically studied in 
the strong coupling limit, and the results are given in the 
Appendix. Due to a frustration induced by the NNN hop- 
ping term, the long-wavelength magnon-mode energy is 
suppressed, and for the simple cubic lattice it vanishes at 
t' = 1/2. At and above this critical value the Neel tem- 
perature Tn vanishes, indicating that the Neel state is 
unstable. The new phase which is stabilized beyond this 
point is a F-AF phase, involving antiferromagnetic order- 
ing of spins in planes (say parallel to the x-y plane), and 
ferromagnetic alignment of spins along the z direction. 
Jill Hence t' = 1/2 marks the phase boundary between 
the AF phase with ordering wavevector Q = (jr, it, it) and 
the F-AF phase with Q = (tt, n, 0). 

For large but finite U, there is a relative stiffening of 
the magnon modes due to the 3rd neighbour ferromag- 
netic coupling induced by the NNN hopping. This re- 
flects a reduction in the degree of frustration, so that a 
slightly higher t' value is required to suppress the magnon 
velocity to zero. From a study of the [/-dependence of 
this critical t' value, the magnetic phase diagram of the 
three dimensional Hubbard model with NNN hopping has 
been obtained recently. jl(| 



V. CONCLUSIONS 



In conclusion, we have studied the spectral function of 
transverse spin fluctuations in an antiferromagnet using 
a unified approach which includes the contributions from 
both single-particle excitations and collective magnon ex- 
citations. For the NN hopping model in two dimensions, 
the integrated spectral weight was studied in the whole U 
range, and the collective magnon contribution was found 
to be dominant for U > 2.5, so that an effective spin 
description of the AF state is appropriate down to a sur- 
prisingly low U value. In the weak coupling limit the 
sharp fall-off of the gapless magnon contribution essen- 
tially leaves only the single-particle contribution having 
a minimum-energy threshold of 2A, so that the spectrum 
of transverse spin fluctuations effectively acquires a gap 
at low energies. The evolution of the spectral function 
with increasing NNN hopping, which reduces the energy 
gap, shows a shift of the spectral function towards lower 
energy due to magnon softening, and also a significant 
rise at low energy due to single-particle excitations. 



G 



APPENDIX 



Hubbard model with NNN hopping — strong coupling limit 

In this section we describe for completeness the AF properties of the half-filled Hubbard model with NNN hopping 
in the strong coupling limit. NNN hopping introduces a frustration in the system through the competing NNN 
AF interaction which enhances spin fluctuations and destabilizes the Neel state. Focussing on the enhancement of 
transverse spin fluctuations at the RPA level by NNN hopping, we examine the spectrum of the collective (magnon) 
excitations, and their contribution to the quantum spin-fluctuation correction <5msF to the sublattice magnetization 
in two dimensions, and the reduction in the Neel temperature Tn in three dimensions. 

As discussed earlier, the quasiparticle amplitudes ai C(T (±) and &k(r(±) which form the eigenvectors of the HF Hamilto- 
nian remain unchanged by the NNN hopping t' , provided the charge gap is finite, which is very much so in the strong 
coupling limit. Therefore the only change in y°{qoj) arises from the change in the quasiparticle energy expression 
given in Eq. (4), which appear in the energy denominators in Eq. (8), and we obtain, for the AA matrix element, for 
example 



[x°(q^)]AA 



E 

k 



E 



V^ 2 



V^ 2 



1 kTe a k-qie 



A 2 

2 



a kte a k-qi 



k-q + ( £ k-q 

2 



A 2 



(11) 



Substituting a 2 .^ = a 2 .^ « 1 — e 2 ./4A 2 and a^ e = a^xe ~ e k/4A 2 in the strong coupling limit, expanding the 
denominator in powers of i/A, t' / A, w/A, and systematically retaining terms only up to order i 2 /A 2 and t 2 /A 2 , 
we obtain after performing the k-sums in two dimensions for a square lattice with Ek e k = 4i 2 , Ek e k — 4i 2 , 
Ek e k<4-q = u ' 2 cos 9x cos q y , 



[x°(q^)]AA = ^ 



i 

u 



At 2 



A 2 



2t 



1- -rs + -Tt(1 

A 2 A 2 



2J 



cosg x cos(j a ) - — 



2t 2 
A 2 " 



(1 — cos q x cosq y ) 



(12) 



where 2A = mil w (1 — 2t 2 /A 2 )U and J = At 2 /U. Similarly evaluating the other matrix elements, [O, we obtain 



[1 - tfx°(qw)] 



A 2 



j' 



-fl - V ) + — 



7q 



7q 
(l-7q) 



2.7 



(13) 



where 7 q = (cos^ + cosq y )/2 and 7 q = cosq x cosq y . Here J = At 2 /U and J' = At 2 /U are the NN and NNN spin 
couplings in the equivalent Heisenberg model, and the NNN term J'(l— 7') directly leads to a softening of the magnon 
mode energies. Substituting in the RPA expression, we finally obtain for the transverse spin propagator 



[x- 



[qu) 



1 (2J 

2 V ufc 



rd-7 q ) 

-7q 



2 J 



-7q 

^(1-7') 



2.7 



1 



CJ — W q + IT] LU + LU q 



11] 



(14) 



where the magnon-mode energy w q is given by 

-*)}'-* (15, 

In the long wavelength limit (g — » 0), with 7 q = 
cos cos <7 y w (1 - g 2 /2)(l - g 2 /2) = 1 - g 2 /2, and 
7 q = (cos q x + cosq y )/2 «1- <7 2 /4, the magnon energy 



reduces to: 

uj cl = V2Jq(l-—) (16) 

showing the strong softening of low-energy modes by the 
NNN hopping. The spin-wave velocity vanishes in the 
limit J'/J — > 1/2. The magnon density of states 
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FIG. 8. The magnon density of states for different values 
of the ratio J'/ J. 

evaluated from Eq. (15) is shown in Fig. 8 for different 
values of the ratio J' /J. NNN hopping clearly softens 
the magnon energies, and transfers the magnon spectral 
weight from the high-energy to the low-energy region. 

The strong softening of the magnon-mode energies sug- 
gests an enhancement in the transverse spin fluctuations. 
To examine this we evaluate the transverse spin cor- 
relations as described earlier. From Eq. (14) for 
the transverse spin propagator, after performing the fre- 
quency integral we obtain the local transverse spin cor- 
relations 
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FIG. 9. The rapid increase in the spin-fluctuation cor- 
rection to sublattice magnetization with the frustrating NNN 
spin coupling J'. 
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The spin-fluctuation correction to sublattice magnetiza- 
tion is then obtained from HT 



FIG. 10. The rapid decrease in the Neel temperature for 
the simple cubic lattice with the frustrating NNN spin cou- 
pling J'. T N /J = 0.989 for J" = 0. 



(5msF = 



(s+s- 



s-s^ 



)RPA 



(s+s- 



- 1 



/RPA 



(18) 



where the denominator (S + S —S S + )bpa is precisely 1 
for the 5 = 1/2 system due to the commutation relation 
[S^S - ] = 2S Z . HH The spin-fluctuation correction to 
sublattice magnetization, evaluated from Eqs. (17), (18) 
is shown in Fig. 9, showing the rapid rise in transverse 
spin fluctuations with the frustrating NNN spin coupling 
J'. 



[x- + (q-)] = ^(- 
2 \u), 



i-^CW,)-^ 

-7q 



We now consider the reduction in the Neel tempera- 
ture in three dimensions due to the frustrating NNN spin 
coupling. For this purpose we consider the NNN hopping 
model on a simple cubic lattice. In this case the lattice 
free-particle energies are 

£k = — 2t(cos k x + cos k y + cos k z ), 

e'k = — 4i'(cos k x cos k y + cos k y cos k z + cos k z cos k x ) , (19) 

and simple extension of the earlier treatment for the two- 
dimensional case leads to the following result for the 
transverse spin propagator at the RPA level 



1 _ ML (l _ ~M + JO. 



"7q 



1 



1 



IT) 



(20) 
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where 7 q = (cosq x + cosq y + cosq z )/3 and 7 q = 
(cos q x cos q y + cosg y cosg z + cosg z cosg;c)/3. The 
magnon-mode energy w q is given by 



uJ 



3J 



^0 \>\ 



-,1/2 



(21) 



For small q, with 7 q » 1 — q 2 /3 and 7 q w 1 — g 2 /6, the 
magnon energy reduces to 



Wo 



Jg ( 1 - 



J 



1/2 



(22) 



which vanishes in the limit J' / J — *■ 1/4 due to the frus- 
tration effect of the NNN coupling J'. The softening of 
the low-energy magnon spectrum has a bearing on the 
Neel temperature, as discussed below. 

Within the renormalized spin-fluctuation theory, |^] 
the Neel temperature Tn is obtained from the isotropy 
condition (S + S~ +S~ S+)t=t n = §5(5+1). For the NN 
coupling model (J' = 0), the Neel temperature was ob- 
tained earlier as Tn = z J ^^-^ fgp for the general case 
of spin S and z nearest neighbors on a hypercubic lat- 
tice. [0 For the simple cubic lattice the spin-fluctuation 
factor / SF = £ q l/(l - 7 q ) = 1-517, and for S = 1/2 
this reduces to Tn/J = 0.989. Extending this analysis 
to the present case, from the expression for the magnon 
propagator in Eq. (20) we obtain 



-In — — r- 



E 

q 



1 - 2 f( 1 -7 q ) 



{i-¥(i-7 q )r 



'q 



(23) 



The Neel temperature, evaluated from the above equa- 
tion, is shown in Fig. 10 as a function of J'. The rapid 
reduction of with J' and the vanishing at J' = J/4 
is due to the enhancement of transverse spin fluctua- 
tions arising from the frustration- induced softening of the 
long- wavelength, low-energy magnon modes. The insta- 
bility at J 1 = J/A is towards a F- AF phase with ordering 



wavevector Q = (7r, tt, 0), involving antiferromagnetic or- 
dering of spins in planes (say parallel to the x-y plane), 
and ferromagnetic alignment along the z direction, fllq] 
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